Superheavy Nuclei: Relativistic Mean Field Outlook. 
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The analysis of quasiparticle spectra in heaviest A ~ 250 nuclei with spectroscopic data provides 
an additional constraint for the choice of effective interaction for the description of superheavy nuclei. 
It strongly suggest that only the parametrizations of the relativistic mean field Lagrangian which 
predict Z = 120 and N = 172 as shell closures are reliable for superheavy nuclei. The influence of the 
central depression in the density distribution of spherical superheavy nuclei on the shell structure is 
studied. Large central depression produces large shell gaps at Z — 120 and N = 172. The shell gaps 
at Z = 126 and N = 184 are favored by a flat density distribution in the central part of nucleus. 
It is shown that approximate particle number projection (PNP) by means of the Lipkin-Nogami 
method removes pairing collapse seen at these gaps in the calculations without PNP. 

PACS numbers: 21.60.Cs, 21.60.Jz, 27.90.+b, 21.10.Pc, 21.10.Ft, 21.10.Gv 



The possible existence of shell-stabilized superheavy 
nuclei, predicted with realistic nuclear potentials and the 
macroscopic-microscopic (MM) method (see references 
quoted in Ref. has been a driving force behind exper- 
imental and theoretical efforts to investigate superheavy 
nuclei. These investigations pose a number of challenges. 
The recent discovery of elements with Z = 115 [|J and 
Z = 116 Q clearly shows great progress on the experi- 
mental side, but also indicates difficulties in the investi- 
gation of nuclei with low production cross sections and 
analyses based on few events. 

The theoretical challenges are also considerable since 
different theoretical methods predict different spherical 
shell closures such as Z = 114, 120, 126 for protons and 
N = 172,184 for neutrons [J]. The largest variation 
in the predictions of shell closures appear in the self- 
consistent calculations based either on non-relativistic 
(Skyrme and Gogny ||) or relativistic (relativistic mean 
field [RMF] @) density functionals. Unfortunately, the 
properties of known superheavy nuclei in their ground 
states do not allow to discriminate between these pre- 
dictions Q- This is due to the fact that known region 
of superheavy elements is dominated by a decay, but a- 
decays occur between neighbouring nuclei and their half- 
lives are only insignificantly modified by shells effects [1| . 

The part of these variations is definetely related to 
the fact that there is a large variety of the parametriza- 
tions for self-consistent models, but for many of them 
even the reliability of describing conventional nuclei is 
poorly known. In addition, self-consistent calculations 
have been confronted with experiment to a lesser degree 
and for a smaller number of physical observables (mainly 
binding energies and quantities related to their deriva- 
tives) as compared with MM method. In particular, 
little attention has been paid to single-particle degrees 
of freedom within the self-consistent models. However, 



the accuracy of predictions of spherical shell closures de- 
pends sensitively on the accuracy of describing the single- 
particle energies, which becomes especially important for 
superheavy nuclei where the level density is high. 

Sect. U shows how the study of quasiparticle states in 
odd- mass nuclei of the A ~ 250 region within the cranked 
relativistic Hartree-Bogoliubov (CRHB) theory con- 
strains the choice of the RMF parametrization for the 
study of superheavy nuclei. The implications of such in- 
vestigation for these nuclei are also discussed (Sect. fl}. 
The influence of central depression in density distribu- 
tion on the shell structure of spherical superheavy nuclei 
is discussed in Sect. [II] The importance of particle num- 
ber projection for the description of pairing properties of 
spherical superheavy nuclei is shown in Sect. Hill Finally, 
Sect. HVl summarizes main conclusions. 
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I. THE A ~ 250 MASS REGION TEST 

The investigation of the single-particle states in the 
odd-mass deformed nuclei of the A ~ 250 mass region 
(the heaviest nuclei for which detailed spectroscopic data 
are available) shed additional light on the reliability of 
the predictions of RMF theory on the energies of spher- 
ical subshells responsible for 'magic' numbers in super- 
heavy nuclei. This is because several deformed single- 
particle states experimentally observed in odd nuclei of 
the A ~ 250 region originate from these subshells (see 
Table 3 in Ref. Q). 

The comparison between experimental and mean field 
single-particle states is less ambiguous in deformed nu- 
clei as compared with spherical ones [l(l HH at least at 
low excitation energies, where vibrational admixtures to 
the wave function are small. This is a result of the sur- 
face vibrations being less collective in deformed nuclei 
than in spherical ones since they are more fragmented 
(ToL H2| . As a consequence, the corrections to the ener- 
gies of quasiparticle states in odd nuclei due to particle- 
vibration coupling are less state-dependent in deformed 
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nuclei. 

A proper description of odd nuclei implies the loss of 
the time-reversal symmetry of the mean-field, which is 
broken by the unpaired nucleon. The BCS approxima- 
tion has to be replaced by the Hartree-(Fock-)Bogoliubov 
method, with time-odd mean fields taken into account. 
The breaking of time-reversal symmetry leads to the loss 
of the double degeneracy (Kramer's degeneracy) of the 
quasiparticle states. This requires the use of the signa- 
ture or simplex basis in numerical calculations, thus dou- 
bling the computing task. Furthermore, the breaking of 
the time-reversal symmetry leads to nucleonic currents, 
which cause nuclear magnetism fl3j| . The CRHB(+LN) 
theory @, 03 takes all these effects into account. 

First ever fully self-consistent description of quasipar- 
ticle states in the framework of the RMF theory was 
presented in Ref. [l[ on the example of 249 > 251 Cf and 
249 Bk nuclei. Fig. Q] shows the CRHB results for 251 Cf. 
Although the same set of quasiparticle states as in ex- 
periment appears, the calculated spectra are less dense. 
This is related to the effective mass (Lorentz mass in 
the notation of Ref. [Hj]) of the nucleons at the Fermi 
surface m*(k F )/m. While the experimental density of 
the quasiparticle levels corresponds to m*(kF)/m close 
to one, the low effective mass m*(fcp)/m ~ 0.66 of the 
RMF theory [l6[ leads to a stretching of the energy scale. 
It has been demonstrated for spherical nuclei that the 
particle- vibration coupling brings the average level den- 
sity in closer agreement with experiment 1 101 ] . Similar 
effect is expected in deformed nuclei. 
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FIG. 1: Experimental and theoretical quasiparticle energies 
of neutron states in 251 Cf. Positive and negative energies 
are used for particle and hole states, respectively. Solid and 
dashed lines are used for positive and negative parity states, 
respectively. The symbols 'NL3' and 'NL1' indicate the RMF 
parametrization. From Ref. [lj 



The calculated energies of a number of states are rather 
close to experiment. On the other hand, the energies 
of some states and their relative positions deviate sub- 
stantially from experiment. For example, only NL1 gives 
the correct ground state v [620] 1/2 in 251 Cf, whereas NL3 
gives the i/[615]9/2 (Fig. [TJ. Detailed analysis shows that 
the discrepancies between experiment and calculations 
can be traced back to energies of spherical subshells from 
which deformed states emerge. This allows us to define 
'empirical shifts' to the energies of spherical subshells (see 
Ref. [l[ for details), which, if incorporated, will correct 
the discrepancies between calculations and experiment 
seen for deformed quasiparticle states. These 'empiri- 
cal shifts' are shown in Fig. [5] as the energy difference 
between self-consistent and corrected energies of specific 
subshells. 
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FIG. 2: Proton single-particle states in a 172 120 nucleus. 
Columns 'NL3' and 'NL1' show the states obtained in the 
RMF calculations at spherical shape with the indicated 
parametrizations. The energy of the I113/2 state in the NL1 
parametrization is set to be equal to that in NL3, which 
means that the energies of all states in NL1 (last column) are 
increased by 0.78 MeV. The columns 'NL3 CO r' and 'NLl cor ' 
show how the spectra are modified if empirical shifts were 
introduced based on discrepancies between calculations and 
experiment for quasiparticle spectra in deformed 249 Bk. Solid 
and dashed lines are used for positive and negative parity 
states. Spherical gaps at Z — 114 and Z — 120 are indicated. 
From Ref. Q 



In the NL1 and NL3 parametrizations, the energies of 
the spherical subshells, from which the deformed states 
in the vicinity of the Fermi level of the A ~ 250 nu- 
clei emerge, are described with an accuracy better than 
0.5 MeV for most of the subshells (see Fig. [2] in the 
present manuscript and Fig. 28 in Ref. [l[: 'empirical 
shifts', i.e. corrections, for single-particle energies are in- 
dicated in both figures). The discrepancies (in the range 
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of 0.6 - 1.0 MeV) are larger for the nlh 9/2 (NL3, NL1), 
^l*n/2 (NL3), ^lji5/2 (NL1) and v>2g$/ 2 (NL3) spherical 
subshells. Considering that the RMF parametrizations 
were fitted only to bulk properties of spherical nuclei this 
level of agreement is good. 

In contrast, the accuracy of the description of single- 
particle states is unsatisfactory in the NLSH and NL-RAl 
parametrizations, where 'empirical shifts' to the energies 
of some spherical subshells are much larger than in NL1 
and NL3. NL-SH and NL-RAl are the only RMF sets 
indicating Z — 114 as a magic proton number p7l.[l8j. In 
the light of present results, these parametrizations should 
be considered as unreliable. 



Single-particle energies [MeV] 
Protons Neutrons 



-8 



2 S U 




120 



V 2 U. 

2 lu- 
H,, 



(a) 



114 



i 

a 
a 



§ 
i 

c 
c 



z 



z 



3s,. 



IK 



2h, 



4s 



184 



3d,, 



3d,, 



172 



(b) 




-l 



-li 



-13 



-15 



FIG. 3: Single-particle spectra of the 292 120 nucleus. 
Columns 'DD-ME2', 'DD-ME1' and 'NL3' show the states 
obtained in the RMF calculations at spherical shape with the 
indicated parametrizations. Fhe column 'NL3 C or' shows the 
spectra modified by empirical shifts. Solid and dashed lines 
are used for positive and negative parity states. 

It is interesting to compare the results obtained 
through non-linear parametrizations of the RMF La- 
grangian (NL1 and NL3) with the ones based on the 
parametrizations which include an explicit density de- 
pendence of the meson-nucleon couplings (DD-ME1 fl9j 
and DD-ME2 [13]). The latter parametrizations provide 
an improved description of asymmetric nuclear matter, 
nucleon matter, and nuclei far from stability [2(|. How- 
ever, the single-particle spectra of the 292 120 nucleus ob- 
tained with DD-ME1 and DD-ME2 are similar to the 
ones seen in NL3 (Fig. [3]) . When compared with the NL3 
spectra corrected for 'empirical shifts' (columns NL3 cor 
in Fig. [3J one can see that the DD-ME2 and DD-ME1 
parametrizations do not remove the problems of the de- 
scription of proton 7rl/ig/2 state (Fig. [3^) and provide 
worse description of neutron vlin/2 state as compared 
with NL3 set (see Ref. [l[ for more detailed discussion of 
these states). 



The measured and calculated energies of the single- 
particle states at normal deformation provide constraints 
on the spherical shell gaps of superheavy nuclei. Such 
analysis restricts the choice of the RMF parametrizations 
only to those which predict Z = 120 and N = 172 as shell 
closures in superheavy nuclei [l[ . In general, since the ac- 
curacy of the description of the v^S\i 2 state is unknown, 
we cannot also exclude the existence of the N = 184 gap 
[l[. Similar analysis to Ref. [1] in non-relativistic models 
would restrict the choice of effective forces. However, it is 
already clear that the SkI4 Skyrme force which predicts 
Z = 114 shell gap can be ruled out since it provides poor 
description of the spin-orbit splittings [16j . Thus, one 
can conclude that non-relativistic theories suggest the ex- 
istence of shell gaps (not necessary doubly shell closures) 
at Z = 120, 126 and N = 172, 184 (see Ref. [H| and 
references quoted therein). The role of self-consistency 
effects in the appearance of these shell gaps is discussed 
in the next section. 



II. SELF-CONSISTENCY EFFECTS IN 
SUPERHEAVY SPHERICAL NUCLEI 

Self-consistent microscopic calculations find a central 
depression in the nuclear density distribution [Tol |22| . 
which generates a wine-bottle nucleonic potential. The 
influence of this depression on the shell structure of 
spherical superheavy nuclei has been studied in Ref. [2lJ 
within the RMF theory without pairing. 

The underlying microscopic mechanism for an appear- 
ance of this central density depression is illustrated in 
Fig.[U The starting point of this consideration is the den- 
sity distribution in 208 Pb, which is nearly flat in central 
region of nucleus. Its charge distribution is well described 
by the RMF theory [2^, [24[ . It is general feature of all 
nuclear structure models that on going from 208 Pb to 
spherical superheavy nuclei in the region around Z = 120 
and N = 172, the ground state configurations are built 
first by the occupation of the group of high-j subshells 
(neutron vliu/ 2 , ^lji5/2, and v2g^/ 2 an d proton 7rli 13 / 2 
and irlh 9 / 2 ) and then by the occupation of the group 
of medium- j subshells (neutron v2g 9 / 2 and proton 2/7/2 
and 2/5/2), see Fig. |U The high-j subshells are local- 
ized mostly near the surface, whereas the low-j subshells 
have a more central localization. As a consequence of this 
grouping of high-j and medium- j subshells above 208 Pb, 
the density is added mostly in the surface region which 
leads to the appearance of the central density depression 
in the nuclei around the Z = 120, N = 172 system, see 
Fig. [4] 

On the contrary, the group of low-j subshells (neutron 
vZd^iii ^3^3/2 and ^4s!/ 2 and proton 7r3p 3 /2 and 7r3p!/ 2 ) 
is filled on going from the Z — 120, N — 172 system to 
the Z — 126, N — 184 system. Since filling up a low— j 
group with nucleons increases the density near the center, 
the density distribution in the central part of nucleus is 
nearly flat in the latter system. 
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The magic gaps of the wine-bottle (the case of cen- 
tral depression in density distribution) and flat-bottom 
(the case of flat density distribution) nucleonic potentials 
are different. This is illustrated in Fig. where starting 
from the ground state configuration of 292 120i72 nucleus, 
having a central depression (configuration 'g-s'), a flat- 
ter density distribution in the central part of nucleus is 
generated by exciting particles from high-j subshells to 
low- J subshells (configuration 'exc-s') [2lJ. While wine- 
bottle potential (conf. 'g-s') is characterized by the large 
Z = 120 and N = 172 large gaps, we see the appear- 
ance of the Z = 126 proton gap and the shrinking of 
the Z = 120 shell gap in the flat-bottom potential (conf. 
'exc-s'). To a lesser extend, the N = 172 neutron gap 
decreases and the N = 184 gap increases in the latter 
potential. 

Due to the isovector force, which tries to keep the neu- 
tron and proton density profiles alike, there is a mutual 
enhancement of the Z = 120 and N — 172 gaps, both 
being favored by the wine-bottle potential, and of the 
Z = 126 and N — 184 gaps, both favored by the flat bot- 
tom potential. For the same reason the gaps are smaller 
for the combination Z = 126 and N — 172, and the 
Z = 120 gap does not develop for the N = 184 systems. 
This behavior is not expected to depend much on the 
density functional chosen (see discussion in Ref. [2lj|). 

Considering that Z = 120, 126 and N = 172, 184 are 
the particle numbers which appear as the candidates for 
'magic' particle numbers in self-consistent theories, it is 
clear that the magnitude of the central depression in den- 
sity distribution is an important factor in defining 'magic' 
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FIG. 4: The evolution of proton and neutron densities with 
the changes of proton and neutron numbers. Arrows indicate 
the group of single-particle subshells which become occupied 
with the change of the nucleon number. The figure is based 
on the results of spherical RMF calculations without pairing 
employing the NL3 [25| parametrization. The shaded back- 
ground is used for nuclei located beyond the proton-drip line. 
If the indicated configuration is not lowest in energ y, i ts exci- 
tation energy (in MeV) is given by E*. From Ref. [21J 



shell gaps in spherical superheavy nuclei. This magnitude 
is correlated with effective mass to* /to of microscopic 
theory: more pronounced central depression develops for 
the density functionals with low effective mass [l6| . This 
feature may be understood as follows [2l| . In the surface 
region, m* /to changes from low value (< 1) in the inte- 
rior to 1 in the exterior. Classically, nucleons with given 
kinetic energy are more likely to be found in regions with 
high effective mass than in the regions with low one be- 
cause they travel with lower speed. This is reflected by 
the Thomas-Fermi expression for the nucleonic density 
p oc [2m*(e F -V)fl 2 . The increase of the effective mass 
in the surface region favors the transfer of mass from the 
center there, which makes the above discussed polariza- 
tion mechanism of the high-j subshells more effective for 
functionals with low effective mass. 

All experimentally known nuclei with Z > 100 are ex- 
pected to be deformed [2(| [13] ■ The deformation leads 
to a more equal distribution of the single-particle states 
emerging from the high-j and low-j spherical subshells 
(see, for example, the Nilsson diagrams in Figs. 3-4 in 
Ref. [2gj |) than for spherical shape. Thus, the de nsity 
profile of a deformed nucleus is relatively flat [U l29j |. 
strongly resembling the one of phenomenological poten- 
tials. In addition, the density profile variations as a func- 
tion of particle number are less drastic than in spherical 
nuclei. These features together with the fact that the 
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FIG. 5: Single-particle spectra of the ground state (indicated 
as 'g-s') and the excited (indicated as 'exc-s') configurations in 
the 292 120i72 system obtained in the RMF calculations with 
the NL3 force. Solid and dashed lines are used for positive and 
negative parity, respectively. Solid and open circles indicate 
the occupied and empty subshells, respectively. In the ground 
state, all subshells below Z = 120 and N = 184 are fully 
occupied. In the excited configuration, only 12 particles are 
excited from the subshell flj'15/2: 4 particles still reside in this 
subshell. The spherical shell gaps of interest are indicated. 
From Ref. [H 
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single-particle energies of the deformed nuclei in heavy 
actinide region have been carefully fitted in the phe- 
nomenological potentials explains the success of the shell 
correction method [30, 31] in the description of known su- 
perheavy nuclei. However, this method neglects the self- 
consistent rearrangement of single-particle levels due to 
the appearance of a central depression in spherical super- 
heavy nuclei. Thus the predictions of the magic numbers 
for superheavy nuclei within the shell correction method 
should be considered with caution. 



III. PAIRING CORRELATIONS IN 
SUPERHEAVY NUCLEI. 

Almost all published self-consistent calculations for 
superheavy nuclei performed either in an approximate 
HF+BCS or full Hartree-Fock-Bogoliubov (HFB) frame- 
works show the collapse of pairing at large shell gaps 
in spherical superheavy nuclei (see, for example, Refs. 
fl7l.l32. 33]. Principal shortcoming of these calculations is 
the fact that neither the BCS nor the HFB wave functions 
are the eigenstates of particle number operator p], |34| . 
The best way to deal with this problem would be to per- 
form an exact particle number projection before the vari- 
ation [341 ] . but this is very time-consuming for realistic 
interactions. In this section, the importance of particle 
number projection on the properties of spherical super- 
heavy nuclei is studied employing the CRHB theory [§] 
with and without approximate particle number projec- 
tion by means of the Lipkin-Nogami (LN) method. 

The CRHB(+LN) calculations are performed using 
the NL3 parametrization for the RMF Lagrangian and 
Gogny D1S force in the pairing channel. The scaling 
factors / of the Gogny D1S force are selected as fol- 
lows: / = 1.0 in the calculations without LN (CRHB) 
and / = 0.864 in the calculations with LN (CRHB+LN). 
These scaling factors provide good description of the mo- 
ments of inertia of rotational bands in the A ~ 250 mass 
region [l[. 

Figs. [Sh.,b compare the calculated pairing energies 
Epairing = — ^Tr(AK) for Z = 120 isotopes obtained in 
the CRHB and CRHB+LN calculations. The Z = 120 
gap is large in the vicinity of N = 172 (see Sect.|H|) which 
leads to the collapse of proton pairing in the CRHB cal- 
culations (Fig. [SJa) . With increasing neutron number, the 
densities in the central part of nuclei become flatter (see 
Sect. [TTJ) leading to shrinking of the Z = 120 gap. Be- 
cause of that the proton pairing shows up at N = 182 
and increases in absolute magnitude with the increase of 
N (Fig. EJb). Because neutron N = 172 and N = 184 
gaps are smaller than the Z = 120 gap (Fig. [3]), the pair- 
ing collapse in neutron subsystem is seen only at these 
neutron numbers (Fig.[SjD). Similar collapse of pairing is 
seen at 'magic' shell closures in the CRHB calculations 
for the N = 172 and N = 184 isotones (not shown here). 
On the contrary, no pairing collapse is observed in the 
CRHB+LN calculations for these chains of nuclei (Figs. 



[6^,,c). In addition, smaller variations of E pa i r i ng as a 
function of particle number are seen in the CRHB+LN 
calculations as compared with the CRHB ones. 

It was suggested in Ref. [33[ to use the fact that the 
pairing energies vanish at closed shells as a fingerprint 
of the shell gaps. Present studies in the CRHB+LN 
framework do not support this suggestion. One should 
also note that the LN method maybe less reliable in the 
regime of weak pairing [35| typical for large shell gaps, 
so more detailed investigation of the pairing in spherical 
superheavy nuclei in the formalism with exact particle 
number projection is needed. 
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FIG. 6: Pairing energies E pa i T i ng as a function of particle 
number obtained in the CRHB calculations with (panels (a) 
and (c)) and without (panel (b)) approximate particle number 
projection by means of the Lipkin-Nogami method for the 
Z — 120 isotopes (panels (a) and (b)) and N = 172 isotones 
(panel (c)). Panel (d) shows the S2 n {Z,N) quantity for the 
TV = 172 isotones obtained in the CRHB calculations with 
and without LN. 

One of the observables used frequently in the search 
of shell gaps in superheavy nuclei is the 5 2n r 2p \(N,Z) 
quantity (called sometimes as two-nucleon shell gap [see 
Ref. [l| for the discussion of this quantity] ) which for the 
neutrons (and similarly for the protons) is defined as 

S 2n (Z, N) = S 2n {Z, N) - S 2n (Z, N + 2) 

This quantity is related to the derivative of the separa- 
tion energy S 2n {Z 1 N), and thus, it is a sensitive indicator 
of the localization of the shell gaps. As follows from the 
comparison of the CRHB and CRHB+LN calculations 
(Figs. [6]: and d), the presence of pairing at shells gaps af- 
fects the S 2n ^ 2p )(N, quantity in the following way: the 
peak in S 2n ^ 2p j(N, Z) becomes more smooth and broad 
and the magnitude of it decreases. 
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IV. CONCLUSION 

The analysis of quasiparticle spectra in heaviest A ~ 
250 nuclei provides an additional constraint for the choice 
of effective interaction for the description of superheavy 
nuclei within specific model. Based on this analysis it 
was concluded that only parametrizations of the RMF 
Lagrangian predicting Z = 120 and N — 184 as shell 
closures provide reasonable description of the spectra of 
odd- mass nuclei. No support for the Z = 114 shell gap 
has been established. One can restrict the predictions of 
self-consistent models (including non-relativistic ones) to 
shell gaps at Z = 120, 126 and N = 172, 184. The inves- 
tigation of self-consistency effects related to central de- 
pression in the density distribution finds important corre- 
lations between these gaps: large central depression pro- 
duces large shell gaps at Z = 120 and N — 172, while 
the shell gaps at Z = 126 and N = 184 are favored 
by a flat density distribution in the central part of nu- 
cleus. However, the magnitude of the central depression 
correlates with an effective mass of nucleon in specific 
model/parametrization. Unfortunately, there is no con- 



sensus on what value an effective mass should have in 
self-consistent theories (see, for example, the discussion 
in Refs. [36j]). Because of existing differences in effec- 
tive mass of nucleon in different models/parametrizations 
(especially pronounced for Skyrme forces), the studies 
of quasiparticle spectra in deformed actinidc nuclei may 
not provide sufficient constraint for localization of magic 
double shell closure (if any) in spherical superheavy nu- 
clei. Approximate particle number projection (PNP) by 
means of the Lipkin-Nogami method removes pairing col- 
lapse seen at large shell gaps in spherical superheavy nu- 
clei in the calculations without PNP. Since closed shell 
nuclei maybe in the regime of weak pairing, more detailed 
studies of the pairing properties of superheavy nuclei in 
the formalism with an exact particle number projection 
are needed. 
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